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Abstract 


The  graph  theory  derived  model  for  the  bonding  topology  in  the  globally 
delocalized  polyhedral  boranes  BgHg2”  and  B^H^-  is  evaluated  by  comparison 
of  the  energies  of  the  core  molecular  orbitals  with  those  obtained  by  the  1962 
LCAO-MO  extended  Huckel  calculations  of  Hoffmann  and  Lipscomb.  Of  particular 
interest  is  how  well  the  complete  graphs  Kg  and  Kj  2  used  in  the  graph  theory 
derived  model  approximate  the  bonding  topologies  of  the  unique  internal  orbitals 
(radial  orbitals)  of  the  octahedron  and  icosahedron,  respectively.  In  the  case 
of  the  BgHg2-  octahedron  the  single  positive  eigenvalue  of  the  Kg  graph 
corresponds  to  the  results  from  the  extended  Huckel  calculations.  In  addition, 
the  parameters  from  the  latter  calculations  indicate  a  ratio  of  0.625  for  the 
overlap  of  the  unique  internal  orbitals  of  the  trans  atom  pairs  relative  to  .£& 
atom  pairs  as  compared  with  unity  implied  by  an  unweighted  Kg  graph.  In  the 
case  of  the  B^H^*”  icosahedron  the  graph  theory  derived  model  is  far  less 
satisfactory  since  the  single  positive  eigenvalue  of  the  K12  graph  disagrees  with 
the  four  bonding  core  molecular  orbitals  (an  Afg  and  three  triply  degenerate 
Tju  molecular  orbitals)  found  by  the  extended  Huckel  calculations  after  removing 
*  the  effect  of  the  mixing  of  core  and  surface  bonding  orbitals  corresponding  to 
the  same  irreducible  representations.  However,  this  core-surface  orbital  mixing 
raises  the  energy  of  the  triply  degenerate  Tju  core  molecular  orbitals  to 
antibonding  levels  so  that  the  graph  theory  derived  model  fortuitously  gives 
correct  skeletal  electron  counts  for  the  regular  icosahedron  despite  this 
fundamental  error. 


1.  INTRODUCTION 


Several  years  ago  we  developed  a  graph-theoretical  approach  for  the  study 
of  the  bonding  topology  in  polyhedral  boranes,  carboranes,  and  metal  clusters.2*3 
Subsequent  work  has  shown  this  approach  to  be  very  effective  in  relating  electron 
count  to  cluster  shape  for  diverse  metal  clusters  using  a  minimum  of  computation. 
Metal  clusters  treated  effectively  by  this  approach  include  post-transition  element 
dusters,4  osmium  carbonyl  dusters,3  gold  clusters, 6*?  platinum  carbonyl 
clusters,6'4  and  rhodium  carbonyl  clusters  having  fused  polyhedra.9'10 

This  graph-theory  derived  method  uses  ideas  originating  from  Huckel  theory.^**2 
However,  in  order  to  eliminate  the  need  for  computation  for  determination  of 
the  signs  of  the  energy  levels,  additional  assumptions  are  introduced  in  order 
to  facilitate  the  estimation  of  energies  at  least  to  the  correct  signs.  This  paper 
examines  in  greater  detail  some  of  these  additional  assumptions  using  energy 
levels  obtained  by  extended  Huckel  calculations  as  a  basis  for  comparison.  Since 
this  comparison  is  facilitated  by  selecting  systems  having  maximum  symmetry 
and  a  minimum  number  of  external  orbitals  to  mix  with  the  polyhedral  core  and 
surface  bonding  orbitals,  the  original  Hoffmann-Lipscomb  LCAO-MO  extended 
Huckel  calculations  on  the  polyhedral  borane  anions  BnHn3  are  used13  selecting 
their  results  for  the  "3N  problem"  of  the  regular  polyhedral  BgHg2-  and  B 1 2 H -j  22-. 

2.  BACKGROUND 

The  topology  of  chemical  bonding  can  be  represented  by  a  graph  in  which 
the  vertices  correspond  to  atoms  or  orbitals  participating  in  the  bonding  and 
the  edges  correspond  to  bonding  relationships.  The  adjacency  matrix  A  of  a 
graph,  such  as  a  graph  representing  chemical  bonding,  can  be  defined  as  follows: 
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i'O  if  H 

1  if  i  and  j  are  connected  by  an  edge  CD 

0  if  i  and  j  are  not  connected  by  an  edge 

The  eigenvalues  of  the  adjacency  matrix  are  obtained  from  the  following 
determinantal  equation: 

1 A  -  xl|  =  0  (2) 


where  I  is  the  unit  matrix  (!;{  =  I  and  l|j  =  0  for  i  4  j). 

The  eigenvalues  of  the  adjacency  matrix  of  the  graph  representing  the  relevant 
chemical  bonding  topology  as  determined  by  equation  2  are  closely  related  to 
the  energy  levels  as  determined  by  Huckei  theory  J  M  4, 15, 16  xhUs  Huckel  theory 
uses  the  secular  equation 

|  H  -  ES(  =  0  .  (3) 

in  which  the  energy  matrix  H  and  overlap  matrix  S  can  be  resolved  into  the  unit 
matrix  I  and  the  adjacency  matrix  A  as  follows: 

H  =  ol  +  BA  (4a) 

S  =  I  +  SA  (4b) 

The  energy  levels  of  the  system  are  related  to  the  eigenvalues  of  the  adjacency 
matrix  A  (equation  2)  as  follows: 


Thus  a  positive  eigenvalue  x  of  A  corresponds  to  a  bonding  orbital  and  a  negative 
eigenvalue  x  corresponds  to  an  antibonding  orbtiai  in  the  corresponding  chemical 
system.  Furthermore,  in  Huckel  theory  S  is  taken  to  be  zero  and  <*  is  used  as 
the  zero  of  energy  (i.e.,  set  to  zero),  in  this  case 

E  =  x6  (6) 

and  the  energy  levels  are  directly  proportional  to  the  eigenvalues  of  the  adjacency 
matrix. 

The  polyhedral  boranes  discussed  in  this  paper  exhibit  globally  delocalized 
bonding.2*3  A  vertex  boron  atom  in  such  systems  uses  three  of  its  four  valence 
orbitals  for  intrapoly hedr a  1  bonding  leaving  one  valence  orbital  as  an  external 
orbital  to  bond  to  the  external  group,  typically  a  monovalent  group  such  as 
hydrogen  or  halogen.  One  of  the  major  triumphs  of  the  graph-theory  derived 
approach  to  the  bonding  topology  in  globally  delocalized  systems  is  the 
demonstration  of  the  close  analogy  between  the  bonding  in  two-dimensional  planar 
aromatic  systems  such  as  benzene  and  that  in  three-dimensional  deltahedral 
boranes  and  carboranes.2  In  this  context  a  deitahedron  is  defined  as  a  polyhedron 
in  which  all  faces  are  triangles.  Thus  consider  a  globally  delocalized  polygonal 
or  deltahedral  system  having  n  vertices.  In  such  a  system  the  three  internal 
orbitals  on  each  vertex  atom  are  divided  into  two  twin  internal  orbitals  (called 
"tangential"  in  some  treatments12)  and  a  unique  internal  orbital  (called  "radial" 
in  some  treatments12).  Pairwise  overlap  between  the  2n  twin  internal  orbitals 
is  responsible  for  the  formation  of  the  polygonal  or  deltahedral  framework  and 
leads  to  the  splitting  of  these  2n  orbitals  Into  n  bonding  and  n  antibonding  orbitals. 


This  bonding  topology  can  be  represented  by  a  disconnected  graph  having  2n 
vertices  corresponding  to  the  2n  twin  internal  orbitals  and  n  Kj  components 
where  a  K2  component  has  two  vertices  with  an  edge  connecting  them.  The 
dimensionality  of  this  bonding  of  the  twin  internal  orbitals  is  one  less  than  the 
dimensionality  of  the  globally  delocalized  system.  Thus  in  the  case  of  the 
two-dimensional  planar  polygonal  systems  such  as  benzene,  the  pairwise  overlap 
of  the  2n  twin  internal  orbitals  leads  to  the  a  -bonding  network  which  may  be 
regarded  as  a  collection  of  n  one-dimensional  bonds  along  the  perimeter  of  the 
polygon  involving  adjacent  pairs  of  polygonal  vertices.  The  n  bonding  and  n 
antibonding  orbitals  correspond  to  the  a  bonding  and  a*  antibonding  orbitals, 
respectively.  In  the  case  of  the  three-dimensional  deltahedral  systems,  the 
pairwise  overlap  of  the  2n  twin  internal  orbitals  results  in  bonding  over  the 
two-dimensional  surface  of  the  deltahedron,  which  may  be  regarded  as 
(topologically)  homeomorphic  to  the  sphere. 

The  equal  numbers  of  bonding  and  antibonding  orbitals  formed  by  pairwise 
overlap  of  the  twin  internal  orbitals  are  supplemented  by  additional  bonding 
and  antibonding  molecular  orbitals  formed  by  global  mutual  overlap  of  the  n 
unique  internal  orbitals.  This  bonding  topology  can  be  represented  by  a  graph 
G  in  which  the  vertices  correspond  to  the  vertex  atoms  of  the  polygon  or 
deltahedron  or  (equivalently)  their  unique  internal  orbitals  and  the  edges  represent 
pairs  of  overlapping  unique  internal  orbitals.  The  relative  energies  of  the 
additional  molecular  orbitals  arising  from  such  overlap  of  the  unique  internal 
orbitals  are  determined  from  the  eigenvalues  x  of  the  adjacency  matrix  A  of 
the  graph  G  (see  equations  2,  5,  and  6,  above).  In  the  case  of  benzene  the  graph 
G  is  the  Cg  graph  (hexagon)  which  has  three  positive  and  three  negative  eigenvalues 
corresponding  to  the  three  *  bonding  and  three  **  antibonding  orbitals, 
respectively.  In  the  case  of  a  globally  delocalized  deltahedron  having  n  vertices 
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such  as  found  in  the  deltahedrai  boranes  BnHn^~  and  the  carboranes  C2Bn_2Hn 
(6  _<_  n  _<_  1 2),  the  graph  C  is  the  complete  graph  Kn  in  which  each  of  the  vertices 
has  an  edge  going  to  every  other  vertex  leading  to  a  total  of  n(n-1)/2  edges.  This 
corresponds  to  an  n-center  bond  at  the  center  (core)  of  the  deitahedron  formed 
by  overlap  of  each  unique  internal  orbital  with  every  other  unique  internal  orbital. 
The  complete  graph  Kn  has  one  positive  eigenvalue,  namely  n-1,  and  n-1  negative 
eigenvalues,  namely  -1  each,  regardless  of  the  value  of  n  indicating  that  the 
n-center  core  bond  in  a  globally  delocalized  deitahedron  leads  to  only  one  new 
bonding  molecular  orbital.  The  sum  of  the  n  bonding  orbitals  arising  from  the 
surface  bonding  of  the  twin  internal  orbitals  and  the  single  bonding  orbital  arising 
from  the  n-center  core  bonding  of  the  unique  internal  orbitals  gives  a  total  of 
n+1  bonding  orbitals  for  globally  delocalized  deltahedra  having  n  vertices.  Filling 
these  n+1  bonding  orbitals  with  electron  pairs  in  the  usual  way  gives  a  total  of 
2n+2  bonding  electrons  in  accord  with  the  observed  number  of  skeletal  electrons 
required  to  form  stable  globally  delocalized  deitahedral  boranes  and  carboranes. 

The  major  objective  of  this  paper  is  to  examine  how  well  the  complete  graph 
Kn  approximates  the  core  bonding  in  globally  delocalized  deltahedrai  systems. 
Thus,  consider  an  octahedral  borane  such  as  BgHg2-.  The  corresponding  complete 
graph  Kg  has  (6)(5)/2  =15  edges.  Among  these  15  edges,  12  represent  overlap 
of  the  unique  internal  orbitals  located  on  adjacent  vertex  atoms  of  the  octahedron, 
namely  a  pair  of  atoms  connected  by  one  of  the  12  edges.  Such  edges  represent 
cis  interactions.  The  remaining  3  edges  of  the  Kg  graph  represent  overlap  of 
the  unique  internal  orbitals  situated  on  one  of  the  three  pairs  of  antipodal  vertices 
of  the  octahedron.  Such  edges  represent  trans  interactions.  Use  of  an  unweighted 
Kg  graph  to  represent  the  core  bonding  topology  in  an  octahedron  gives  equal 
weights  to  the  cis  and  trans  interactions  despite  their  obvious  geometric 
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difference.  This  problem  was  already  recognized  by  Schmidtke  in  1967  in  one 
of  the  early  papers  pointing  out  the  relationship  between  graph  theory  and  Huckel 
theory.^  In  group  theoretical  terms  the  graph  theory  derived  model  uses  the 
symmetric  group  Sg  (i.e.,  the  automorphism  group  of  the  complete  graph  Kg) 
rather  than  its  subgroup  0^  (i.e.,  the  point  group  of  the  octahedron)  to  represent 
the  symmetry  of  the  skeletal  bonding  manifold  of  octahedral  boranes  such  as 
BgHg2~,  which  consists  of  the  6K2  surface  bonding  and  the  Kg  core  bonding. 
The  actual  Of,  rather  than  Sg  symmetry  of  these  systems  will  result  in  partial 
or  complete  removal  of  the  five-fold  degeneracy  of  the  core  antibonding  orbitals 
and  the  six-fold  degeneracies  of  the  surface  bonding  and  antibonding  orbitals. 

A  similar  situation  occurs  in  icosahedrai  boranes  such  as  The 

corresponding  complete  graph  K 12  has  (12)(11)/2  =  66  edges.  Among  these  66 
edges,  30  edges  represent  overlap  of  the  unique  internal  orbitals  located  on 
adjacent  vertex  atoms  of  the  icosahedron,  namely  a  pair  of  atoms  connected 
by  one  of  the  30  edges  and  corresponding  to  the  ortho  interaction  using  carborane 
nomenclature.  An  additional  30  edges  of  Kj2  represent  overlap  of  the  unique 
internal  orbitals  located  on  non-adjacent,  non-antipodal  vertex  atoms  of  the 
icosahedron,  namely  the  meta  interaction.  The  remaining  6  edges  of  Kj2  represent 
overlap  of  the  unique  internal  orbitals  located  on  the  6  pairs  of  antipodal  vertices 
of  the  icosahedron,  namely  the  para  interaction.  Use  of  an  unweighted  K|2  graph 
to  represent  the  core  bonding  topology  of  an  icosahedron  gives  equai  weights 
to  the  ortho,  meta.  and  para  relationships  despite  their  obvious  geometric 
difference.  The  graph  theory  derived  model  thus  uses  the  symmetric  group  S12 
of  order  12!  =  479,001,600  rather  than  the  icosahedrai  group  1^  of  order  120  to 
represent  the  symmetry  of  the  skeletal  bonding  manifold  of  icosahedrai  boranes 
such  as  Bi2Hi2^”'  which  consists  of  the  I2K2  surface  bonding  and  the  Kj2  core 


bonding. 

Reduction  of  the  symmetry  from  that  of  the  symmetric  groups  Sn  (n  =  6  or 
12,  respectively)  to  the  polyhedral  groups  and  1^  for  octahedral  and  icosahedral 
boranes,  respectively,  not  only  splits  some  molecular  orbitals  of  high  degeneracies 
but  also  leads  to  some  mixing  of  the  core  and  surface  bonding.  Thus  in  the  case 
of  octahedral  boranes  under  0^  symmetry  the  core  bonding  has  the  representation 
A1g  +  Tlu  +  Eg  and  the  surface  bonding  has  the  representation  Tju  +  T2g  +  T2U 
+  Tig.  The  Tty  core  an<*  surface  bonding  molecular  orbitals  can  therefore  overlap 
leading  to  corresponding  lowering  and  raising  of  the  energies  of  the  two  molecular 
orbitals  (designated  as  AE(Tju)).  Therefore,  the  Tju  orbital  energies  determined 
by  computational  methods  based  on  0^  symmetry,  including  the  extended  Huckel 
method,  do  not  correspond  to  pure  core  or  surface  bonding  but  also  include  the 
interaction  AE(T]U),  which  represents  another  variable  which  needs  to  be  evaluated 
before  the  graph  theory  derived  methods  can  be  compared  with  various 
computations.  Thus  in  an  octahedral  borane  such  as  BgHg2-  under  Oh  symmetry, 
only  the  Ajg  and  Eg  molecular  orbitals  can  represent  pure  core  bonding.  In 
octanedral  metal  clusters  in  which  external  bonding  orbitals  clearly  play  a  major 
role  (e.g.,  CogfCOJi^-  studied  by  Mingos^},  even  the  Ajg  and  Eg  molecular 
orbitals  do  not  represent  pure  core  bonding  because  of  interaction  with  other 
orbitals  of  Ajg  and  Eg  symmetries  arising  from  the  external  bonding. 

A  similar  situation  occurs  for  icosahedral  boranes.  Under  lh  symmetry  their 
core  bonding  has  the  representation  Ag  +  Tju  +  T2U  +  Hg  and  their  surface  bonding 
has  the  representation  Tju  +  Hg  +  Gu  +  Gg  +  Hu  +  T|g.  The  Tju  and  Hg  core 
and  surface  bonding  orbitals  mix  with  energies  AE(Tju)  and  AE(Hg),  respectively, 
leaving  only  Ag  and  T2u  as  pure  core  bonding  orbitals. 

This  paper  compares  the  graph-theory  derived  and  Huckel  theoretical  models 
for  BfiHg2-  and  Bj2Hj2^~  by  determining  which  edge  weights  for  the  corresponding 


Kn  graphs  retaining  the  polyhedral  symmetries  Of,  and  1^,  respectively,  reproduce 
the  energies  for  the  pure  core  orbitals  (A^g  and  Eg  for  BgHg2-  and  Ag  and  T2U 
for  B-j2h122 “)  calculated  by  the  extended  Huckel  method13  avoiding  externa! 
orbitals  (i.e.,  the  "3N"  set  of  Hoffmann  and  Lipscomb13)  in  order  to  eliminate 
core-external  orbital  mixing  effects.  These  edge  weights  are  then  used  to 
calculate  hypothetical  energies  for  the  impure  core  orbitals  (Tju  for  BgHg2** 
and  Tju  and  Hg  for  B^H^2”1  ,n  the  absence  of  mixing  (unperturbed  core  energies). 
Comparison  of  these  hypothetical  energies  with  the  actual  energies  from  the 
extended  Huckel  calculations13  allows  evaluation  of  the  core-surface  mixing 
energies  (A.E(Tiu)  for  BgHg2**  and  AE(Tiu)  and  A  E(Hg)  for  B^H^2-)*  The 
assumptions  of  the  simple  graph  theory  derived  bonding  model  can  then  be 
evaluated  using  the  following  observations: 

(1)  The  single  positive  eigenvalue  of  the  complete  Kn  graphs  should  correspond 
to  the  completely  symmetric  (Ajg  or  Ag)  core  orbital  being  the  only  bonding 
orbital.  All  other  core  orbitals  should  have  antibonding  unperturbed  core  energies 
corresponding  to  the  n-1  negative  Kn  eigenvalues. 

(2)  The  weights  of  the  edges  in  the  polyhedral !y  weighted  Kn  graphs  representing 

non-adjacent  vertices  in  the  octahedron  or  icosahedron  should  be  close  to  the 
unit  weights  taken  for  the  edges  in  the  Kn  graphs  representing  adjacent  vertices 
in  the  polyhedron.  - 

The  first  observation  is  the  important  observation  affecting  applicability  of  the 
graph  theory  derived  methods.  Deviations  of  edge  weights  from  unity  are  not 
serious  if  they  do  not  create  more  than  one  positive  eigenvalue  for  the  core 
bonding. 

3.  THE  OCTAHEDRON 

Consider  an  octahedrally  weighted  Kg  complete  graph  having  12  edges  of 
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unit  weight  corresponding  to  the  octahedron  edges  and  the  remaining  3  edges 
of  weight  t  corresponding  to  the  three  octahedron  antipodal  pairs.  The  spectrum 
of  this  graph  can  be  determined  by  symmetry  factoring?0  using  a  thr  *e-foid 
axis  (Figure  1).  The  symmetric  branch  (A  in  Figwe  1,  top)  gives  the  eigenvalue 
4+t  corresponding  to  the  Ajg  molecular  orbital  as  well  as  one  of  the  -t  eigenvalues 
of  the  triply  degenerate  Tju  molecular  orbital.  The  doubly  degenerate  E  branch 
gives  the  other  two  -t  eigenvalues  of  the  triply  degenerate  T]u  molecular  orbital 
as  well  as  the  two  -2+t  eigenvalues  of  the  doubly  degenerate  Eg  molecular  orbital. 
Note  that  any  positive  value  of  t  is  sufficient  to  lead  to  only  one  positive 
eigenvalue,  namely  the  4+t  eigenvalue  of  the  Ajg  orbital,  and  five  negative 
eigenvalues,  namely  the  -t  eigenvalues  of  the  triply  degenerate  Tju  orbitals 
and  the  -2+t  eigenvalues  of  the  doubly  degenerate  Eg  orbital.  This  indicates 
that  any  positive  trans  interaction  in  an  octahedron  gives  the  same  distribution 
of  bonding  (1)  and  antibonding  (5)  orbitals  as  the  unweighted  Kg  graph  used  in 
the  graph  theory  derived  model?  for  BgHg?",  C2B4H5,  and  globally  delocalized 
octahedral  metal  clusters.  Thus  for  octahedral  boranes  and  metal  clusters  the 
success  of  this  model  is  not  very  dependent  upon  the  value  taken  for  t.  Note 
also  that  setting  t  =  0  gives  the  spectrum  of  the  octahedron  (+4,  0,  0,  0,  -2,  -2) 
whereas  setting  t  =  1  gives  the  spectrum  of  the  Kg  complete  graph  (+5,  -t,  -1, 

The  octahedron  has  two  pure  core  bonding  molecular  orbitals,  namely  Ajg 
and  Egr  in  the  absence  of  mixing  with  external  orbitals  (see  above).  The  computed 
values  for  the  energies  of  these  two  orbitals  substituted  into  equation  6  generate 
two  equations  in  t  and  S;  these  equations  can  be  solved  to  determine  these 
parameters  provided  that  the  zero  energy  level  separating  bonding  and  antibonding 
orbitals  is  known.  Table  1  summarizes  the  results  of  such  calculations  using 
the  energies  computed  by  Hoffmann  and  Lipscomb1 3  for  the  "3N  orbitals"  of 
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BgHg2".  Thus  the  computed  energies  for  the  Ajg  and  Eg  core  bonding  molecular 
orbitals  generate  the  following  equations  in  t  and  S: 


A 1  gJ 

(4+06  =  2.969 

(7a) 

EgJ 

(-2+06  =  -0.884 

(7b) 

Solving  these  two  equations  gives  t  =  0.62S  and  6  =  0.642.  These  values  of  t 
and  6  can  then  be  used  to  calculate  a  hypothetical  value  for  the  Tju  core  molecular 
orbital  in  the  absence  of  core-surface  interaction  using  the  following  equation: 

Tlu:  (-0.625X0.642)  =  -0.401  (8) 

This  compares  with  an  energy  of  -0.829  found  in  the  Hoffmanrt-Llpscomb 
calculation  for  the  Tju  core  orbital  thereby  indicating  the  interaction  parameter 
4  E(T |  u)  to  be  -0.401  -  (-0.829)  *  0.428.  Applying  this  correction  to  the  energy 
of  the  Hoffmann-Lipscomb  calculation  of  the  T surface  orbital  gives  energies 
of  0.595  and  0.493  for  the  bonding  surface  orbitals  (T^u  +  T2g)  and  -0.416  and 
-0.671  for  the  antibonding  surface  orbitals  (T2U  +  Tjg)«  This  distribution  of  the 
energies  of  the  surface  bonding  orbitals  is  crudely  consistent  with  the  prediction 
by  the  graph  theory  derived  method*  of  six  bonding  and  six  antibonding  surface 
orbitals  with  the  same  absolute  values  for  the  energies  arising  from  the  spectra 
of  six  equivalent  disconnected  K  2  graphs. 

Attempts  have  been  made  with  much  less  success  to  apply  this  method  to 
more  recent  computations  on  octahedral  boranes  and  metal  clusters.  Armstrong, 
Perkins,  and  Stewart*^  report  self-consistent  molecular  orbital  calculations 
on  highly  symmetrical  borane  cage  anions  including  BgHg*".  The  zero  energy 


level  of  their  calculation  is  not  defined  so  thata  as  well  as  B  and  t  are  unknowns; 
thus  equation  5  must  be  used  rather  than  equation  6.  This  leads  to  an 
underdetermined  system  with  only  two  equations  to  determine  three  unknowns; 
solution  of  this  system  to  obtain  meaningful  numbers  is  not  feasible. 

The  best  example  of  a  relevant  calculation  on  a  globally  delocalized  octahedral 
metal  cluster  is  the  Wolfsberg-Helmholtz  molecular  orbital  calculation  by  Mingos*9 
on  Cos(CO)i44~.  However,  because  of  the  large  number  of  external  orbitals 
involved  in  this  and  other  metal  clusters,  even  the  Ajg  and  Eg  molecular  orbitals 
do  not  represent  pure  core  bonding.  Thus  the  CogfCO)^^  system  has  three 
molecular  orbitals  of  Ajg  symmetry  and  four  molecular  orbitals  of  Eg  symmetry 
thereby  precluding  an  analysis  of  the  type  given  above  for  the  Hoffmann-Lipscomb 
BgHg2-  computation. 

«.  THE  ICOSAHEDRON 

Consider  an  icosahedrally  weighted  K  J2  complete  graph  having  30  edges  of 
unit  weight  corresponding  to  the  icosahedron  edges,  30  edges  of  weight  m 
corresponding  to  the  meta  interaction  of  non-adjacent  non-antipodal  vertex  pairs, 
and  6  edges  of  weight  p  corresponding  to  the  para  interaction  of  the  6  icosahedron 
antipodal  pairs.  The  spectrum  of  this  graph  can  be  determined  by  symmetry 
factoring  using  first  a  threerfold  axis  and  then  a  two-fold  axis  (Figure  t,  bottom). 
The  fully  symmetric  branch  (AC)  gives  the  eigenvalue  5+5m+p  corresponding 
to  the  Ag  molecular  orbital  as  well  as  one  of  the  -1-m+p  eigenvalues  of  the 
quintuply  degenerate  Hg  molecular  orbitals.  The  remaining  four  -1-m+p 
eigenvalues  of  the  quintuply  degenerate  Hg  molecular  orbitals  arise  from  the 
two  isolated  vertices  of  the  doubly  degenerate  EC  branch.  The  AU  branch  and 
the  doubly  degenerate  EU  branch  generate  a  quadratic  equation  whose  roots 
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give  the  eigenvalues  of  the  triply  degenerate  Tju  and  T2U  molecular  orbitals. 
Note  that  the  Tju  molecular  orbital  has  a  positive  eigenvalue  unless  p  >/S  (1-m). 
Thus  with  most  likely  values  of  the  edge  weights  m  and  p  the  icosahedrally 
weighted  K 1 2  graph  (Figure  1)  has  four  positive  eigenvalues  (the  Ag  and  triply 
degenerate  Tju  orbitals)  rather  than  only  the  single  positive  eigenvalue 
characteristic  of  the  unweighted  Kj2  graph.  Note  also  that  setting  m=0  and  p=0 
gives  the  spectrum  of  the  icosahedron  (+5,  +  /5,  +/5,  +  / 5,  -1,  -1,  -1,  -1,  -1,  -/5, 
-A,  -A)  whereas  setting  m=1  and  p=1  gives  the  spectrum  of  the  K 12  complete 
graph  (+11,  -1,  -1,  -1,  -1,  -1.  -1.  -1,  -1,  -1,  -1,  -1). 

The  icosahedron  like  the  octahedron  has  two  pure  core  bonding  molecular 
orbitals,  namely  Ag  and  T2u.  The  computed  values  for  their  energies  substituted 
in  equation  6  generate  two  equations  in  the  three  parameters  m,  p,  and  0.  The 
system  is  therefore  underdetermined  and  an  additional  relationship  between 
m,  p,  and  8  must  be  used  in  order  to  obtain  unique  values  for  the  three  parameters. 
Instead  of  assuming  such  a  relationship,  consider  the  minimum  value  of  p  to  be 
zero  and  the  maximum  value  of  p  to  be  m.  Then  plot  against  p  in  this  range 
the  corresponding  values  of  the  parameters  m  and  8  (Figure  2)  and  the  molecular 
orbital  energy  levels  adjusted  by  removal  of  the  core-surface  interaction  (Figure 
3)  using  the  data  from  the  Hoffmann-Lipscomb  3N  extended  Huckel  computations 
on  Bj2Hi2^-*  The  values  of  0  and  the  core  molecular  orbital  energies  are  seen 
from  Figures  2  and  3  to  be  rather  insensitive  to  the  actual  value  of  p  taken. 

By  far  the  most  important  conclusion  from  these  data  is  that  for  all  likely 
values  of  p  the  core  bonding  in  icosahedral  B12H122”  leads  to  four  positive 
eigenvalues  (Ajg  +  T-ju)  rather  than  the  single  positive  eigenvalue  obtained  by 
using  the  K12  complete  graph  for  the  core  bonding  topology.  The  graph  theory 
derived  model  *  is  thus  incorrect  for  the  icosahedron  since  it  leads  to  the  wrong 
number  of  positive  eigenvalues  for  the  core  bonding,  namely  one  rather  than 


four.  Nevertheless,  the  graph  theory  derived  model  works  for  icosahedrai  boranes 
and  carboranes  because  this  fundamental  error  is  corrected  by  mixing  the  bonding 
Tju  core  orbitals  with  the  bonding  Tju  surface  orbitals.  Thus  the  mixing  energy 
AEtTiii)  is  so  large  (1.1  to  1.6  on  the  scale  of  Table  2)  that  the  core  Tju  orbitals 
are  raised  to  antibonding  energies  by  the  core-surface  interaction  so  that  only 
one  core  bonding  orbital,  namely  the  Ag,  orbital,  remains  after  this  interaction. 
This  is  why  the  graph  theory  derived  model  for  icosahedrai  loads  to 

the  experimentally  observed  skeletal  electron  count  even  though  the  assumption 
of  the  K  complete  graph  for  the  core  bonding  topology  leads  to  the  incorrect 
number  of  core  bonding  orbitals  (before  core-surface  mixing). 

Important  experimental  support  for  these  ideas  is  found  in  the  centered 
icosahedrai  gold  cluster  Aui3Cl2(PMe2Ph)io3+  (ref.  22)  in  which  there  is  no 
surface  bonding  comparable  to  that  of  Bj2Hl2^~  because  the  p  orbitals  of  the 
12  peripheral  gold  atoms  in  this  cluster  have  energies  too  high  to  function  as 
twin  internal  orbitals  for  surface  bonding.*'7  Turning  off  the  surface  bonding 
in  this  way  in  the  centered  icosahedrai  Auj3Cl2(PMe2Ph)io3+  leads  to  four  core 
bonding  orbitals  in  accord  with  the  energy  levels  of  the  icosahedrally  distorted 
K 12  graph  for  any  reasonable  levels  of  the  parameters  m  and  p  (Figure  1  bottom 
and  Figure  3).  This  point  may  also  relate  to  the  general  observation*  that  the 
overlap  topology  of  the  n-1  peripheral  gold  atoms  in  a  centered  Aun  cluster  leads 
to  electron  counts  corresponding  not  to  those  expected  for  a  complete 

graph  but  instead  to  that  of  the  polyhedron  formed  by  the  peripheral  gold  atoms. 

5.  SUMMARY 

Determination  of  the  eigenvalues  of  the  octahedrally  weighted  Kg  complete 
graph  indicates  one  positive  eigenvalue  for  any  positive  weighting  (t)  of  the 
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antipodal  (trans)  interactions  relative  to  unit  adjacent  fcis)  interactions.  This 
indicates  that  the  unweighted  Kg  graph  used  in  the  graph  theory  derived  model 
for  the  octahedral  boranes,  carboranes,  and  metal  clusters2  is  adequate  for 
determining  the  correct  number  of  core  bonding  orbitals  and  therefore  the 
corresponding  electron  counts.  The  Hoffmann-Lipscomb  calculations1 3  on  BfiHfi2" 
correspond  to  a  value  of  0.625  for  this  parameter  t. 

A  similar  determination  of  the  eigenvalues  of  the  icosahedratly  weighted 
K72  complete  graph  indicates  four  positive  eigenvalues  for  any  likely  values 
for  the  weighting  parameters  m  and  p.  This  is  in  disagreement  with  the  graph 
theory  derived  model 2  for  icosahedral  boranes  and  carboranes  which  uses  the 
unweighted  K|2  graph  for  the  core  bonding  leading  to  only  a  single  positive 
eigenvalue  corresponding  to  a  single  core  bonding  orbital.  However,  core-surface 
bonding  interactions  raise  the  energies  of  three  of  the  four  core  bonding  orbitals 
to  antibonding  levels  so  that  despite  this  major  inaccuracy  the  graph  theory  derived 
model  fortuitously  gives  correct  electron  counts  for  icosahedral  boranes. 

The  analysis  in  this  paper  shows  that  the  graph  theory  derived  model2  describes 
a  bonding  topology  which  is  accurate  enough  to  lead  to  correct  electron  counts 
in  octahedral  clusters  but  leads  to  difficulties  in  icosahedral  clusters.  However, 
most  globally  delocalized  metal  clusters  fortunately  are  octahedra  to  which 
this  simple  bonding  approach  can  be  reliably  applied.  It  thus  appears  that  the 
graph  theory  derived  models  will  continue  to  be  useful  in  understanding  the 
relationship  between  shape  and  electron  count  in  the  large  variety  of  known 
polyhedral  boranes,  carboranes,  and  metal  clusters. 
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TABLE  1 


MOLECULAR  ORBITAL  ENERGY  PARAMETERS 
FOR  OCTAHEDRAL  B6H62~ 


Molecular 

Orbital 

-Cflrg-Bflndlna 

Tlu 

Eg 

Surface  Bonding 
Tlu 

T2g  . 

t2u 

T1g 

Other  Parameters 

6 

t 

AE(T1u) 


Calculated 
by  Hoffmann 
and,  Lipscomb3  (3N) 

2.969 

-0.829 

-0.884 

1.023 

0.493 

-0.416 

-0.671 

0.642 

0.625 

0.428 


Adjusted  by 
Removal  of  Core- 
Surface  Interaction 

2.969 

-0.401 

-0.884 

0.595 

0.493 

-0.416 

-0.671 

0.642 

0.625 

0 


a)  From  R.  Hoffmann  and  W.N.  Lipscomb,  J.  Chem.  Phys..  JLSL, 
2179  (1962). 


0.024  0.164  0.322 


Figure  1:  Symmetry  factoring  schemes  for  the  octahedrally  weighted  Kg  graph 
(top)  and  the  isocahedrally  weighted  K  ]  2  graph  (bottom).  The 
designations  of  the  branches  correspond  to  that  used  in  R.B.  King. 
Theor.  Chim.  Acta.  223  (1977).  The  symmetry  factoring  of  the 
octahedrally  weighted  Kg  graph  uses  a  three-fold  axis  (C3)  leading 
to  an  A  branch  and  a  doubly  degenerate  E  branch.  The  symmetry 
factoring  of  the  icosahedrally  weighted  K  \ 2  graph  first  uses  a  three-fold 
axis  (C3)  similarly  leading  to  A  and  doubly  degenerate  E  branches, 
which  in  turn  use  two-fold  axes  (C2)  for  further  symmetry  factoring 
into  symmetric  (C)  and  antisymmetric  (U)  branches;  the  resulting  four 
branches  are  designated  as  AC,  AU,  EC,  and  EU. 

Figure  2:  Plot  of  the  parameters  m  and  8  versus  p  for  the  Hoffmann-Lipscomb 

3N  extended  Huckel  computations  on  Bi2Hi22-» 

« 

Figure  3;  Plot  of  the  energies  of  the  core  molecular  orbitals  Ag  +  T]u  +  T2U  +  Hg 
versus  p  in  the  absence  of  core-surface  interactions  for  the  Tju  and 
Hg  orbitals  for  the  Hoffmann-Lipscomb  3N  extended  Huckel  computations 
on  Bi2Hj2^~* 
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